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Abstract 



Circular orbits of spinning test particles and their stabihty in Schwarzschild-like backgrounds are in- 
OO ' vestigated. For these space-times the equations of motion admit solutions representing circular orbits with 

particles spins being constant and normal to the plane of orbits. For the de Sitter background the orbits are 
always stable with particle velocity and momentum being co-linear along them. The world-line deviation 
fj ' equations for particles of the same spin-to-mass ratios are solved and the resulting deviation vectors are 

^^, used to study the stability of orbits. It is shown that the orbits are stable against radial perturbations. The 

general criterion for stability against normal perturbations is obtained. Explicit calculations are performed 
kjq' in the case of the Schwarzschild space-time leading to the conclusion that the orbits are stable, 
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2 ■ 1 Introduction 

rf-\ ■ Orbits of particles in curved space-times have long been a topic of interest in physics and astronomy. In 

J ' I relativistic theories of motion, orbits of particles around astronomical objects have been studied extensively in 

^— s . the past and have resulted in a plethora of phenomena ranging from the perihelion precession and bending of 

("^ I light to more recent ones like the so-called clock effect [T]. Nearly circular orbits have had a especial role in 

this regard. Their application in theories dealing with accretion discs [2] (see also [3] and references therein) 
is an example showing their importance. Some of the objects moving on such orbits may have some internal 
angular momentum. In trying to model the motion of such rotating objects treated as test particles, one is lead 
to the Mathisson-Papapetrou-Dixon (MPD) equations [11[51[S1[7]- According to these equations, in general the 
t^ . trajectories of a spinning particle moving in a curved background differ from the geodesies of the background 

space-time as a result of a coupling between the particle spin and the space-time curvature and furthermore, 
its momentum and velocity are not in general parallel. 

Being a system of highly coupled differential equations, the MPD equations are usually difficult to solve 
and only a few analytical solutions to them are known, e.g. the one introduced in |H] in which the motion in the 
equatorial plane of a Kerr black hole with a fixed particle spin was investigated. Following the ideas presented 
in the latter reference, solutions describing circular or nearly circular orbits in various black hole space-times 
have been found. Different aspects of such orbits have been studied in the case of Lense-Thirring space-time 
in [9] , in the case of Schwarzschild space-time in [lOl [TTJ [T2j [131 HI] (also in [15] for non-circular orbits) , in the 
case of Kerr black hole in [H [HI [^ [H [20l [211 [H [23l [24] , in the case of Reissner-Nordstrom black hole in 
[25] , and in the case of Weyl space-time in [55] . An alternative framework has been used in [57] [55] and [SS] 
to study orbits of spinning particles in Schwarzschild and Kerr-Newman space-times respectively. 

An interesting question is whether these orbits are stable. For the case of (spinless) test particles this 
question has been considered in [301 131] where deviation of circular geodesies in static spherically symmetric 
space-times was studied. The aim of the present work is to investigate this (in)stability for orbits of spinning 
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particles. The technique we use in this regard is to study the growth of perturbations to a given orbit. 
These perturbations are essentiaUy deviations from a fiducial orbit to adjacent ones. Unbounded growth of 
a perturbation signals instability of the orbit. A recent application of this technique may be found in j32) 
in which the geodesic deviation equation is used to study the stability of circular orbits of test particles 
in Einstein-Gauss-Bonnet gravity. In fact as it has been argued in |33| the deviation vectors may be used 
to define Lyapunov indices in general relativity for either geodesic or non-geodesies flows. The method of 
Lyapunov indices has been applied to study the stability of circular or non-circular orbits of spinning particles 
in [21 1101 1111 115| . In particular in |11| it was shown that using different supplementary conditions may result in 
different stability behaviour of spinning particles orbiting the Schwarzchild black hole. The problem of chaotic 
motion of spinning particles in Schwarzschild space-time has also been studied in j34j but in that work the 
deviation equations were not used and a numerical scheme based on a time series method was used instead. 
Chaotic motion of spinning particles around the Schwarzschild black hole has also been studied in [35 from a 
different perspective. The case of Kerr balck hole has been considered in [55] . 

To study perturbations of spinning particles orbits, the geodesic deviation equation is no longer the appro- 
priate equation, instead we should use its generalization as presented in [37] or [38] . A similar set of equations 
have been obtained in [39] for the case of charged spinning particles. In this work we deploy the machinery 
developed in [3S] to obtain the deviations. 

As it has been argued in [3D] the test particle approach leading to equations like the MPD equations, is 
valid only when the particles spins are small compared to the characteristic length of the background they are 
moving in. Thus throughout this work we assume the particles specific spins (i.e., their spin-to-mass ratios) to 
be small compared with the mass of the black hole. 

The paper is organized as follows. In the next section we give a brief review of the MPD equations and 
show how a solution describing a circular orbit can be found in a rather general Schwarzschild-like space-time. 
Although circular orbits in particular spherically symmetric backgrounds have been studied in detail in some 
of the above mentioned references, that section serves both as a rather unified manner review of the main 
results essential for the subsequent sections and to fix the notation. The latter is of importance because a 
variety of notations have been used in the literature. These include the way in which the gauge is fixed and 
the particle spin tensor is related to the spin vector. The case of de Sitter space-time is simpler and is studied 
in a separate section as a warm up. Then we solve the world-line deviation equations for circular orbits and 
use them to study the behaviour of the perturbations. Even though our focus of attention is on orbits in 
Schwarzschild space-time, we obtain some general criteria under which the orbits are (un) stable in a general 
Schwarzschild-like space-time. The paper is concluded with a discussion of the results. 

2 Circular Orbits 

The basic equations governing the motion of a spinning test particle are [7] 

^ = -Ir'u.pV^s-^ (1) 



^^ p'^v'^-p^v^, (2) 

P^.s^" = (3) 

in which v^ — ^ , ,P^, s^'^ are the particle velocity, momentum, and spin respectively, -^ — w"Vq with V^ 
being the covariant derivative, and 

represents the space-time curvature tensor. The particle spin is related to its spin vector by [S] 

^ £^"^"^5/3 (4) 



■myf—g 



with e^'^'a^ being the alternating symbol. In this work we adopt e*''^* = -1-1. It should be noted that a different 
defining equation is sometimes used in this regard (see e.g. equation Al of [34]). As a consequence of these 



equations one can show that the particle mass and spin are constant 

p^iP^ = —m= const., (5) 

Sfj,s^ ~ s^ ~ const. (6) 

We choose the following gauge fixing relation [7] 

Pt^v^ = -m (7) 

in which the particle instantaneous zero- momentum and zero- velocity frames are simultaneous 41; . This is the 
gauge in which the velocity-momentum relation [S] |41) 



l2 



'"'' -'-P'+ Jf!^tZ'"\..s P' (8) 



which is itself a consequence of the equations of motion, takes the simplest form. 

We consider the motion of a spinning test particle in a Schwarzschild-like background space-time described 
in the coordinate system {t, r, 0, (p) by the line element 

ds^ ^ -f{r)dt'^ + -^ dr^ + 7-2 {de^ + sin^ Od^^ ) (9) 

in which the explicit form of /(r) is subject to the Einstein field equations. This metric corresponds to 
Schwarzschild space-time for /(r) = 1 — 2M^ to Reissner-Nordstrom space-time for /(r) = 1 — ^^ + ^, to 
de Sitter space-time /(r) = 1 + ^7-^, to Schwarzschild-de Sitter space-time for /(r) = 1 — ^^ + ^r^, and to 
Reissner-Nordstrom- de Sitter space-time for /(r) = 1 — 2M _(_ ^ _(_ A ^2^ 

Inspired by the results of 0, we consider a spinning particle moving on a circular orbit r = a in the 
equatorial plane 6 — ^ with a velocity v^ = (w*, 0,0, u"^). We take the particle spin vector perpendicular to 
this plane. Thus we set 

s^ = 1(0,0,5,0) (10) 

r 

and look for a consistent solution 

p^(r) = (p*(r),p'-(T),0,/(r)) (11) 

of the equations of motion ([T]) and ^. Note that Eqs. ((TU| and (|lll) are consistent with Eq. ^ by construction. 
Now making use of Eq. ([4]), we can insert this ansatz into Eq. ([T]), to obtain 

'^P J t r , ^J <h r |^r,\ 

— vp H w P , (12 

^ = /(^-^z;V + ai.V + ^(-/'W/-F/'A')), (13) 

'€ - -^V~g.y, (14) 



where / = /(a), /' - '^^^'^ 



dr 



s „ 

and s = — . Similarly, equation ^ leads to 
m 



a dr 2a 



v'p'- + v'l'p\ (15) 



§d^ - -^v'p' + sv^p^ + v^p'-v'pt (16) 

fa dr 2a 



sa^ = -sz^V + wV- (17) 

dr 



Solving Eqs. dH]) , dH]) , ([12]) , and ^7^ for the components of p^ lead to 

p* = const., (18) 

p^ = 0, (19) 

p"* = const. (20) 

Now the remaining two equations, P^ and (ITC| . together with Eqs. ([5]) and ([7]) reduce to the following system 
of algebraic equations 

/p*'-aV' = m^ (21) 

/wV-a^wV = n^, (22) 

y w p - aw^"^ = — w> ^^' P^, (23) 

sf 

— vV-swV = wV-w*/ (24) 

2a 

from which we can obtain the non- vanishing components of v^ and p'^. 

3 Motion in de Sitter space-time 

To proceed further, we first consider the simplest case where u'^ and p'^ are parallel, i.e. p^ = mv^ but s 7^ 0. 
Taking this into account, the above system of algebraic equations are consistent only for an /(r) of the form 
of /(r) = 1 + Ar^, A being a constant. This corresponds to the well-known de Sitter space-time. Thus we set 

f{r)^l + K\ (25) 

This restriction means that there is no solution describing a circular orbit with velocity and momentum being 
parallel in other Schwarzschild-like space-times. From the fact that de Sitter space-time is homogeneous we 
already expected w^ and p^ to be parallel. This is easily confirmed by the above equations resulting in 

p^ = TOi'^ = m I 1,0,0,±W- I . (26) 

The angular velocity of particle is independent of the orbit it moves on. 

We now aim to investigate the stability of these orbits. To this end, we consider the behaviour of pertur- 
bations in different directions. Unbounded growth of a perturbation is a signature of instability. A suitable 
framework to study stabilities is the world-line deviation equations developed in [38]. However, in this case the 
momentum and velocity four- vectors are co-linear and the equations of motion reduce to the geodesic equation. 
In fact for maximally symmetric space-times we have 

Rfj,iyal3 ~ C{gfj,ag,^f3 — gii/^g^a) (27) 

where C is constant. Inserting this into Eq. ([T|) returns 

Dp^ 



Dt 



-Cw"s^,. 



On the other hand by combining the above relation with Eq. ^ and making use of Eq. ([3]) we arrive at the 
previously mentioned result p'^ — mv'^. Thus 

Dpf" C 

Ut m 

i.e. particles move along geodesies and their momentum and spin are parallel. Hence we can deploy the 
well-known geodesic deviation equation instead of the world-line deviation equations. It reads 

-^-^ = ^R^o..pn''v"v^ (28) 



where n^ = ^jr- is the deviation between adjacent trajectories x^,x^ + Xn^ . Applying this equation to the 



above circular orbits results in 



% = -//V^ + 2a/.^^-.n^ (30) 

dr'^ ar ar 

= 0, (31) 



dr"^ 



l!L ^ „f„0^ (32) 

dr^ a dr 

where 

fc = ^(/'' + //")(t'*)'-K' + /)(«^)'- 

By integrating the first, the third, and the last equation in the above system of equations we obtain 

dn^ f 

^ = -^v'n^ + C, (33) 

n^ = Csi + Cs, (34) 

^ = --«*«'■ + C4 (35) 

dr a 

in which Ci,C2,C3,C4 are constants. The second equation above shows that the circular orbits considered 
here are stable against normal (i.e., in the ^-direction) perturbations. To see the effect of radial perturbations, 
we first insert the values of /, w* , u"^ into the expression for k to get k = ^. We then insert this back into Eq. 
(pO| and solve it for n'' with the aid of Eqs. (|33|) and ([34]) to obtain 

n'' ^ N cos{2v'''t) + I (36) 

in which N, I are constants. This leads us to the conclusion that the orbits are stable against radial perturba- 
tions. 

4 Motion in other Schwarzschild-like space-times 

Relaxing the simplifying assumption of the previous section, we now seek general solutions p'^ , v^ to Eqs. (1211) - 
(|24p. In general these equations admit two sets of solutions which, roughly speaking, correspond to two senses of 
rotation with respect to a given spin direction. Labelling the solutions by (p^,p!^, w^, w^j!.) and (p^,p_, v*_^v'_) 
respectively, we have 

A^^e^, (37) 



^ 2(,sV-^ - 2){s^f + sa^f" - iasf ± ah) 
^^ ay^{s^af" + s^f' ~4:a±sh) 



-sf + saf"±h t 



2(2a-s2/') 
2{2a~Pn 



where 



2/(2a - s2/')P* - a^i-sf + saf" ± h)p'l' ' 



h = \l -if'^P - Qaf'Pf" + Pa^f"^ + 8f'a + 2f'^s-^f", 



v\ (39) 

(40) 



b± 



f'2-2 



-2s2aV + 8afs^f" - 8f'a + 2/ 
-16/sV" + 4/,sV"' + 16/ ± 2s(/' - a/")/^ 



2/'^sV" 



An interesting property of these solutions is that under s — > — s the absolute values of v'^ and p'^ change in 
general. This is due to the spin-dependent character of the force exerted on the particles. Also note that 
the + and — solutions do not represent physically distinct situations. This can be easily seen from the above 
expressions in which s ^^ —s interchanges -f and — solutions (up to an overall minus sign). In fact only two 
of the four possible situations {s = ±^ sense of rotation = clokwise, counter — clockwise} are independent 
because there is no preferred up or down direction in a spherically symmetric space. 

To get more physical insight into the nature of the above solutions we consider the case of Schwarzschild 
space-time. The functions v'^ /v* are depicted in Fig. [T]in which the range of s is chosen as \s\ < 0.1 to conform 
with the requirement s <^ M and the range of a is taken as a > 3 since for 2 < a < 3, w* becomes imaginary, 
i.e. there is no circular orbit for this range. The case a = 3 needs special care and is discussed later. 




Figure 1: The angular velocity of particles -4- (left) and -;- (right) in terms of a and s. 



The asymmetry s <-> — s is obvious in these graphs. They also show that for large values of a the two 
solutions coincides. To ease comparison, the absolute value of \v't/v[_\ is depicted in Fig. [2l It coincides with 
the left-hand side plot in Fig. [T]upon taking s — >■ —s. 

We have also depicted the functions w±/w± — p'^/p^. in Fig. [3] It shows that for larger values of a the 
difference between p'^ /p* and v'^ /v* diminishes. This is because a larger r = a results in a value of / closer to 
unity and hence a smaller curvature which in turn results in a smaller spin-curvature coupling. 

To investigate the special case of a = 3 mentioned earlier, we first consider the case of a particle without 
spin. For s = 0, Eq. ([M)) is satisfied identically, and the remaining Eqs. (PT|) - (P5)) result in 



2/ - a/' 



(41) 



which for the Schwarzschild space-time implies a > 3. When s ^ 0, Eqs. pT|l - ((24l) do admit a solution given 

by 



~27s±l 
18(27 -s2)' 



(42) 



^, _ v/6l(27' 



P± 



.s2)(2s3 - 135s ± 3/)(-4s3 + 27s ± 3Z) 
2s(27-s2)(9s2=Fs?-|-486) ' 

^6s(27 - s2)(2s3 - 135s ± 3Z) 
' 18s(27-s2) 



(43) 
(44) 




Figure 2: The function |-^| in terms of a and s. 





Figure 3: The non-co-hnearity -j- r- (left) and -j- r (right) in terms of a and s. 



where I — %/— 24s^ + 1053s^ + 8748 and ± in Eqs. (j42|) and (|44p are chosen according to the sign of s. 

Two further remarks are in order. First, an ansatz corresponding to particles moving on circular orbits 
with the spin vector making an arbitrary angle with the plane of the orbit, say s^ ex (0, s, 0, 0), is not in general 
compatible with the equations of motion. Second, since the Schwarzschild-like space-times are spherically 
symmetric there is no preferred equatorial plane and hence similar circular orbits on other planes are also 
allowed. 



5 Stability of orbits 

In general, orbits of spinning particles are not geodesies of the background they are moving in and we can not 
use the geodesic deviation equation. In this work we use the world-line deviation equations [3 8) 



L)t z z 



Dt 



Dn'^ 



2^' ""'^7 



^a/3 



Dt 



(45) 
(46) 



Sfj.iyJ 



-JnuP" 



(47) 



where j'^ — n°'\/aP'^ smd J^^ = n°'\I ^s^^ ■ The dimensions of various quantities used in these equations are 
shown in Table [T] 



Table 1: The dimensions of various quantities used in world-line deviation equations. Here L abbreviates 
length, T = t,r, and iJj 



Quantity 


s 


p' 


p^ 


n^' 


n* 


f 


,^ 


jrr 


jT^ 


ji'i' 


Dimension 


L 


L 


L^ 


L 


L'> 


L 


L^ 


L' 


L 


L'> 



The above equations are the generalization of geodesic deviation equation to world-line deviations of spin- 
ning particles of the same specific spins. In fact setting s'^'^ = in Eq. (1461) immediately results in 



Dnt" 

^d7 



in which 5 is a constant. Inserting this back into Eq. (HSl) we arrive at Eq. (pS]) . the geodesic deviation 
equation. It is interesting to note that, according to Eq. psj) . even if we neglect the direct effect of spin, i.e. 
the first three terms in the right hand side of that equation, the particles spin still affect the dynamics through 
p^ in the last term which in turn depends on spin. The following additional equations 



PmJ 



= 0, 



Pm 



J"" = - 



Jt^i 



(48) 
(49) 



resulting from the constants of motion Eqs. ([5]) and ([5]) respectively are also useful in this regard. 

We now consider a family of circular orbits and take a typical circular orbit of definite radius in the 
equatorial plane along which the particle spin is constant and orthogonal the plane of orbit, as the fiducial 
path. Deviations from this orbit can be obtained by solving Eqs. (|^5|) and (I46p for n^ but before embarking 
on that, it would be useful to reduce the number of unknowns by using Eqs. (l47l) - (|49)) . From Eq. (|48l) we have 



J : 



Q2p0 



(50) 



and from Eq. (|49l) we have 



Similarly Eq. (|47l) result in 



jr^ ^_ii_jtr^ 



pt 



J'"- = - 



sfp' .t 



ap''- 



J^f 



af 



fp' 
op.p't 



jte 



Now, taking Eqs. ([5n)) - ((M)) into account, from Eq. (H5|) with ^ ~ t,r,d, and (j) we obtain 



dT 



saf'p'l' drf f 



2/ 



,^+^iv^~sv^) 



fp' 



ap^n 



4, J 



(51) 

(52) 
(53) 
(54) 

(55) 



d£_ 

dT 



saff'p't' dn* sff'p* dn'l' saff"'v*p^ 



dT 



-(w* - sv't') 



P_^,2 



2 dT 

- Vf"K 



2a V'^ 



(56) 



df _ fff\, , /(1-/)„VV*' 



^V + ^vY + ' 'I, " v^p' ) n\ (57) 



dr \ 2a a?p' 

2a a 



^^K-s«*)(^^^ -ap^'^*--/j, (58) 

respectively. Similarly from Eq. (1461) with fiv = ir, i0, i0, rO, and 0(/) we obtain 

^^-f'7--/.'+p'^, (60, 

dT dr dr ap^pf 

+ (2/p*w'^ - a/yi;* + asfp'f'v'l' - sff"v^p^)ff (61) 

^ = - {-^v^ - ^) J- - f ^.y + ^^^^ " ^^.V) n^, (62) 

dr V 2 ap"^ y V 2 a J 

fP'''''^'-v^r^+.^f-p^'^ (63) 



a^p't' dr a dr 

respectively. Note that setting nv — rcf) lead to Eq. ((5^ again. 
Eqs. d^, dSni), 0, dSi), and dH]) resuh in 

n' = const. 

which is simply separation between two adjacent circular orbits, and 

j = const., j = const., —— = const., — — = const. 

dr dr 

Thus radial perturbations do not lead to instability. 

The solutions to Eqs. dST]), §0^), §2^, and §^ are given by 

^^(t) = Di sin(u£T) + Z?2 cos(w^r) + D^ exp(K±r) + D4 exp(— k±t) (64) 

and similar expressions for j^{t), J*^ (t), J^^{t). Here Di, D2, -D3, D4 are constants, 

y/a±u*p* 
2am 

a± = (2/ - a/')(«V/3± " 2a2//3±7± + sa//'/3± + 2.s/(l - /)7±), 
-{s^f'^+sa^f" ^3asf'±ah) 



P± 



a{Paf" + s^f' ~4:a±sh) 



7± 



-sf + saf" ± h 



Note that under s -H- — s we have a+ <-> a_. Thus for a± < we have osciUatory sin and cos terms only 
and the orbits are stable. For a± > there are both oscillatory and exponential terms and the orbits are not 
in general stable. Depending on the values of space-time parameters and the particles specific spins s, both 
cases may occur. For the case of Schwarzschild space-time we have depicted a± in terms of a and s in Fig. 21 
According to this figure both solutions corresponds to stable orbits in Schwarzschild space-time, at least for 
the range of parameters chosen here. 





Figure 4: a+ (left) and a_ (right) in terms of a, s. 



6 Discussion 

We have studied circular orbits of spinning particles with their spins perpendicular to the plane of motion 
in spherically symmetric space- times. We have shown that circular orbits of spinning particles with their 
velocity and momentum being parallel are only allowed in de Sitter case. For other cases the particle velocity 
and momentum are not parallel. We have obtained general expressions for the particle angular velocities and 
momenta in terms of the radii of orbits and specific spins. Depending on the sense of rotation with respect to 
the spin direction, there are in general two physically equivalent sets of such expressions for any given radius. 
According to these expressions, for a fixed spin direction and orbit, the absolute values of angular velocities of 
particles orbiting in opposite senses are different. This is due to the spin-dependent force exerted on particles. 
For the case of Schwarzschild space-time, there is no circular orbit of radius smaller than three (measured in 
units of the black hole mass). This is in agreement with the results of pii;. The effects due to spin decrease 
when the radius of the orbits increase because the curvature and hence the spin-curvature coupling decreases 
with radius. 

We have also solved the world-line deviation equations of spinning particles of the same specific spins 
moving on orbits described above. The resulting deviation vectors have been used to show that the orbits 
are stable against radial perturbations. We have determined the general condition under which the orbits 
arc stable against normal perturbations. It shows that in general particles orbiting in opposite senses have 
different deviations which may result in different stability properties. By explicit calculations we have shown 
that the orbits are stable in Schwarzschild space-time. An advantage of our approach is that some rather 
general results can be extracted from the equations without resort to detailed numerical calculations. These 
results might be of interest in studies related to (in)stability of non-geodesic flows. The deviations n^ obtained 
here enable us to compute the corresponding Lyapunov indices. Although the special case studied here, i.e. 
exactly circular orbits with a perpendicular spin may not be realized in real astrophysical situations, it can 
be regarded as a toy model which can be extended to more realistic situations. Possible interesting extensions 
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include a study of stability of spinning particle orbits in a Kerr-type background and charged spinning particles 
in a Reissner-Nordstroni space-time. 
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